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Abstract. We consider a spin-orbit coupled system of particles in an external 
trap that is represented by a deformed harmonic oscillator potential. The spin- 
orbit interaction is a Rashba interaction that does not commute with the trapping 
potential and requires a full numerical treatment in order to obtain the spectrum. 
The effect of a Zeeman term is also considered. Our results demonstrate that 
variable spectral gaps occur as a function of strength of the Rashba interaction 
and deformation of the harmonic trapping potential. This implies that the many- 
body physics of spin-orbit coupled systems can be manipulated by variation of 
the external trapping potentials. 
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1. Introduction 

The last decade is associated with breakthroughs in ultracold atomic and state-of-the- 
art optical lattices experiments, which provided not only extremely important data 
for understanding of the fundamentals of quantum mechanics, but also introduce a 
sophisticated set of tools for investigation of exceptionally pure and tunable quantum 
systems [2 HI 13 H]. Recently, this toolbox has been expanded to also include the 
possibility of applying controllable gauge fields to both bosonic O [6l [71 [3 [9] and 
fermionic systems [13 HI] (see references [T^] or [13] for short reviews). Spin-orbit 
coupling is a prime example of a non-abelian potential that plays an important role 
throughout physics. Examples are the spin-orbit splittings in atomic and nuclear 
spectra, and the distinct effect it imposes on the band structure of solid-state systems 
that has lead to great recent advances in the exploration of materials with robust 
metallic surface states, the so-called topological insulators [Ml [15]. 

In many experiment with ultracold atoms, the external trapping potentials, while 
present, are often quite shallow and can be ignored for many purposes or one can 
include them using a local density approximation. However, some recent experiments 
have demonstrated that tight trapping is not only possible but also a very exciting 
possibility as this brings the physics one can study close to what is known from atomic 
or nuclear structure. Some recent theoretical papers have shown that the trapping 
potential can play an important role in spin-orbit coupled systems [13 HZl 113 dS] ■ 
The cited works consider the case where the trapping potential is isotropic in all three 
dimensions ^] or isotopic in the plane where the spin-orbit coupling acts [T71 (13 dS] . 
In the current presentation we extend this discussion by considering also the case 
where the trapping potential can be deformed. 

Our findings show that deformation of the external potential has a decisive 
effect on the spectral density and can lead to the opening and (near) closing of the 
energy gaps in the single-particle spectrum. In the strongly deformed effectively one- 
dimensional limit the level spacing is completely determined by the shallow trap 
harmonic frequency at low energy. However, for small deformation one can find 
regimes where the spectral gaps will almost close and produce a (quasi)-continuum. 
This implies that deformation can be used as a control parameter for the level 
density. In the case of a many-body system with non-abelian gauge potentials, the 
spectral density plays a decisive role in trapped systems when exploring many-body 
phenomena such as exotic pairing and crossover f23 HI [22l 123 ES], superfluidity 
and condensation [23 HZl |23 123 130], ferromagnetism [3T] or quantum Hall states 
[33 133 1311 133 133 133 133 132] • Our results represent an initial step in exploring these 
interesting questions for spin-orbit coupled systems in deformed traps. 

2. Formalism 

We consider the single-particle problem in a deformed three-dimensional (3D) 
harmonic trap including a Rashba [30] spin-orbit interaction and an applied magnetic 
field. Here we use the Rashba form of the spin-orbit interaction but the formalism 
is completely general and applies to all types of spin-orbit coupling including the 
case where the Rashba-type and Dresselhaus-type [41] contributions are equal (which 
amounts to a term like cr^py in the notation introduced below) which is used in many 
current cold atomic gas experiments studying spin-orbit effects. The Hamiltonian for 
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a particle with mass m experiencing a Rashba and Zeeman term is 

2 1 ^ 2 1 ^2 



TT Px , ^ 2 2 , Py , ^ 2 2 , Pz , ^ 22, / \ Ty f-\\ 

2m 2'^^"^ 2m 2^^y^ 2m 2^^^ +aKKpy-ayP^)-/x-B, (1) 



where are Gy are 2x2 Pauli matrices and is the strength of the Rashba spin-orbit 
coupling. We consider a particle with two internal degrees of freedom, which we label 
as spin up |t) and spin down \X). We choose the effective magnetic field as B = (0, 0, B) 
and wo can write the Zeeman term as —fiBa^ with ctz the remaining Pauli matrix and 
fj, is the effective magnetic moment. For the harmonic trapping potentials we assume 
different trapping frequencies in all directions, uj^, t^y and w^. However, since we 
take the Rashba term to act only in the xy-plane the z-direction effectively decouples 
from the problem. Note that in the case of B = 0, the Hamiltonian is symmetric 
under time-reversal and for a single particle Kramer's theorem dictates a two-fold 
degeneracy. 

The Hamiltonian can be written in the matrix form 

anipy - ipx) Hox + Hoy + f^B J l^Vi / V^; / ' 

where we have introduced the short-hand Hqx = ^ + ^mw^x^ and similarly for Hoy. 
This constitutes an effective two-dimensional (2D) problem. The oscillator energy in 
z-direction can be considered as a parameter that can be included in the; total energy, 
E. The natural basis in which to expand our wave- functions are eigenfunctions of the 
2D harmonic oscillator. Thus, 



V'T = XI «i("x,«3/) I n^,ny,t) (3) 

V'i = X bi{i^x,ny) \n^,ny,i), (4) 

where 

with 63. = and = — ^ . Introducing the standard ladder operators „ and 



al y that obey 



X = 



h , rnhLOx, t X 



^^(4 + ay), Pv = \J ^—{al-ay), (6) 



the Hamiltonian now reads 



h{N. + l) + iNy + l)--^ ^{iV^ial - ay) - ^7(4 - (^pA ^ /Vr 

[^al - ay) + V7(4 - «x)) 7(^x + 5) + (^x + i) + if ) UJ ' U 

(7) 

where we have introduced the notation = a^Ux and Nx = a^ax, the dimensionless 
Rashba coupling /3 = c(r^/ 2fitr^ ^^'^ ^ ~ '^^^ ratio of the harmonic oscillator 
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Figure 1. Energy as function of tlie dimensionless spin-orbit coupling parameter 
13 for tlie case of equal frequencies uj = ujx = with no Zeeman shift (loft panel) 
and including a Zeeman shift of magnitude fiB = hui. 



frequencies is 7 = We will be using hujy as the unit of energy. The linear system 
of equations for the coefficients a and b becomes 



ilii^x + ^)+ny + ^)- e)a{n^, Uy) + /? 



i^nyb{nx,ny - 1) - i^/r^^~+Tb{nx, Uy + I) 



\/l\/nih{nx - 1,%) + y/l\/n.^ + lb[nx + l,ny) 



= 



ili^x + ]j)+ny + ]^)- e)b{n^, Uy) + (3 



i^/n:^a{nx,ny - 1) - iy^Uy + la{nx, Uy + f ) 



+ \/l\/n^a{nx - 1, riy) - ^\Jn^ + \a{n^ + 1, riy) 



= 0, 



(8) 



where eo{nr^,ny) = j{nx + |) + {uy + |). This set of equations cannot be solved 
analytically and one has to resort to numerical methods. In the symmetric case where 
7 = 1, one could also have used a basis based on the solution of the harmonic oscillator 
potential in cylindrical coordinates [T7| [181 [19] . However, this is not appropriate in our 
case since the cylindrical symmetry is broken in the plane for 7^1 and the problem 
is therefore better handled using the Cartesian basis expansion presented above. 



3. Results 



We now present the results of our study of the single-particle spectral structure in 
the presence of a Rashba spin-orbit term and with an external trap that can be 
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deformed. All the results presented here are in the regime < /3 < 1.5. It is a 
simple matter to go to even higher Rashba couplings but it requires the use of a 
bigger single-particle basis in order for all states to properly converge. The results we 
present below have all been obtained by using a basis with about 700 single-particle 
states. However, due to deformation the maximum number of quanta, max{nx) and 
max(nj,), is generally not equal to each other since we cut the basis in energy space, 
i.e. ma.x{eo (nxjUy)) = 7 (max(na;)-|-l/2)-|-(max(ny)-|-l/2) is our cut-off that restricts 
the values of and riy . 

3.1. Symmetric Trapping Potential 

In the left panel in figure [T] we show the single-particle energy spectrum in a two- 
dimensional spherical harmonic oscillator as function of the coupling strength, /3. The 
oscillator degeneracies that are well-known for /3 = are lifted as (5 increases and the 
oscillator shells become more and more mixed. Each level is still doubly degenerate 
due to time reversal symmetry mentioned above. The lowest levels decrease and for 
sufficiently large strengths the behavior approaches a parabolic dependence on /3 that 
can be found in a semi-classical treatment |16j . The many levels are spread out over 
energies and overall they cover rather densely the energy space. In many regions the 
spectrum resembles a continuous distribution where the individual levels cannot be 
distinguished. The deviations from the regular picture seen in the top right-hand 
corner on both panels in figure [l] is due to the numerics and the use of a finite basis 
set. We have numerically checked that this can be remedied by systematic expansion 
of the basis size. 

The key features are found in the low-density pockets of the energy spectrum 
where no levels are found. We see that we have a spectrum that can open and 
close 'super' gaps as a function of the Rashba coupling strength. These regions are 
significant since we would expect that a system with this level density would tend 
to be more stable with respect to the application of weak interactions between the 
particles or perturbations from the background since the larger gaps hinder strong 
mixing of states. In the opposite case where the level density is high we expect to see 
strong mixing and small perturbations could have large effects. What is particularly 
remarkable here is the fact that the Rashba term can almost cancel the effects of the 
trap, i.e. the spacing of trap levels at /3 = of one huj is strongly reduced at /3 ~ 0.2 
for total energies above ~ 5fia;. Indeed these sorts of changes are reflected on the 
properties of a many-body system in a symmetric trap as discussed for the case of 
condensates in references [T7], [TS] and [TS]. 

On the right panel in figure [l] we show again the case of equal frequencies, 
but this time with the inclusion of a Zeeman field of magnitude 
fiB = hu!. The Zeeman field will not influence the oscillator levels but only displace 
the spin components. Note that the absolute ground state of the spectrum now starts 
at zero energy and then decreases. We have opted to keep the same vertical scale on 
both panels in figure [l] in order to make a comparison of the overall spectral density 
at higher energies and therefore the Zeeman shifted ground state level cannot be 
seen. What we find is that the Zeeman split will tend to counteract the effect of the 
Rashba term at small /3 so that we now find a spectrum with large gaps at /3 ^ 0.2. 
However, since the spin splitting sends states up and down in the spectrum it tends to 
generate a sort of 'two-gap' structure around f3 ^ 0.2 with one large and one smaller 
gap, intertwined by a number of densely spaced states. The conclusion is that both 
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Figure 2. Energy as function of dimensionless spin-orbit coupling parameter 
/3 for the case where the oscillator potential is deformed. The left panel has 
7 = — = 2 and the right panel has 7 = 3. 



Rashba and Zeeman can be used as an experimental handle on the density of states in 
the single-particle spectrum as the density and structure can be changed by varying 
one while keeping the other fixed. 

3.2. Deformed Trapping Potentials 

We now consider the case of an oscillator potential that is non-spherical corresponding 
to different frequencies in the x and y directions. In this situation, high degeneracy is 
found for special frequency ratios of deformed harmonic oscillators. These special 
configurations occur when the frequency ratios equal ratios of small integers as 
2/1, 3/2, etc. This is due to the harmonic oscillator spectrum being linear in 
both the frequencies and quantum numbers in the different spatial directions. The 
harmonic oscillator has especially high degeneracy compared to deformations of 
(almost) all other radial shapes. This means that gaps are more likely to occur in 
the spectrum. The implication for pure oscillators with no Rashba or Zeeman terms 
is that degeneracies are periodically occurring as function of the frequency ratio, and 
that degeneracies are washed out in between the highly degenerate points. In figure [2] 
we show the spectrum for frequency ratios 7 = 2 (left panel) and 7 = 8 (right panel) as 
a function of Rashba coupling strength, (3, with no Zeeman field. For these ratios we 
clearly see that the spectrum is going towards a more degenerate form with stronger 
gaps for increasing 7 in comparison to the 7 = 1 case shown in figure [T] Indeed for 
7 = 3 we have to go up to energies of IbHujy and above before we can see the same 
behavior as for lower energies in the left panel of figure [T] When including a Zeeman 
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Figure 3. Same as figure|2]but with ratios 7 = 1.57 in the left panel and 7 = 2.17 
in the right panel. 

term for these ratios (not shown here) we see the same trends as seen from left to 
right in figure [T] i.e. the spectrum becomes slightly more dense overall and there is 
competition between Rashba and Zeeman terms. 

A completely different scenario is displayed in the case where 7 is not a ratio 
of small integers. In figure [3] we show the single-particle spectrum for 7 — 1.57 
(left panel) and 7 = 2.17 (right panel). For these ratios the states are more evenly 
distributed and for 7 = 1.57 we see an almost constant spectral density for energies 
of bhujy and above, while for 7 = 2.17 this is not seen until about lOhujy and above. 
Comparing the left and right-hand panels in figure |3j we see an overall tendency for 
the larger 7 to have a smaller overall density of levels since this is closer to the one- 
dimensional limit that we will return to momentarily. The overall quadratic decrease 
with f3 is still seen as in figures [l] and [2j Comparing the results presented in figure [2] 
for the ratios 7 = 2 and 7 = 3 to those of figure [3] with 7 = 1.57 and 7 — 2.17 
we thus conclude that the deformation of the trap is another experimental way to 
change the spectral structure and density. However, for the latter ratios the influence 
of the Rashba coupling is diminished somewhat as the density changes in a much 
more smooth manner as compared to figures [T] and [2j This implies that the choice 
of deformation is very important when studying the effects of spin-orbit coupling on 
trapped systems. 

3.3. One- dimensional Limit 

Increasing the frequency ratio towards infinity separates the Hamiltonian in a low- 
energy one-dimensional part very weakly coupled to the high-energy one-dimensional 
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Figure 4. Same as figure [2] and [s] for 7 = 5 (left panel) and 7 = 10 (right 
panel). These results approach the limit of an effectively one-dimensional system, 
i.e. 7 > 1. 



part in the other direction. With our conventions, the hmit 7 <C 1 corresponds to a 
strongly confined motion along the a;-direction and a shallow confinement along the y- 
direction. These decoupled one-dimensional equations can be solved analytically with 
the corresponding Rashba couplings. The Pauli matrix for either or ay implies that 
one should use a basis of eigenfunctions for the operator with the small frequency, i.e. 
Wy, since the other direction has a very large oscillator frequency that effectively 
freezes the motion. This linear combination of equal amplitude spin-up and spin- 
down spinors immediately decouple the equations of motion into separate oscillator 
eigenvalue problems. The eigenvectors of are 

*±-L\V (9) 



\±1 

and these can now be used as the two-component basis instead of the spin-up and 
spin-down spinors. The one-dimensional Schrodinger equation then becomes 

^ + -mojy±aRPy-EU±f±{y) = 0, (10) 

where f±{y) is the radial wave function. We rewrite these equations as 



1 \2 1 22 rnai, 
— {py±maR) +^mujyy 



^±f±{y)^o, (11) 



from which the harmonic oscillator solutions can be directly inferred to be 

Eriy ^ -^ma% + tmjyiny + ^) (12) 
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The first term is simply the quadratic decrease of the energy with /3 that we already 
noted above and that can be obtained by semi-classical means [TB]. We can compare 
this dispersion to the one-dimensional solution with no external trapping potential in 
the 2/-direction which is simply 

Ek^ = ^±aRhky, (13) 

where ky is the one-dimensional momentum along the y-direction. These two 
dispersion relations look very different but can be reconciled by recalling that the 
virial theorem tells us that (p^) cx mhuiyriy where py — pyzL maR and we thus recover 
a momentum-dependence of the energy with quadratic and linear terms. 

In figure |4] we show the spectral structure as one approaches the one-dimensional 
limit. In the left panel of figure |4] we have 7 = 5 and in the right panel 7 = 10. The 
low-energy eigenvalues reduce to the equidistant harmonic oscillator spectrum with 
the same frequency but shifted quadratically with the spin-orbit coupling strength. 
The remarkably simple emerging feature is that the lowest part of the spectrum for 
an even modest deformation already resembles the equidistant one-dimcnsional limit. 
This can be clearly seen by making a comparison of figure |4] to figures [2] and [3] As 
the deformation increases an increasing part of the low-energy spectrum approaches 
the one-dimensional limit. This feature can be understood from the weak coupling of 
two oscillators with very different frequencies. The perturbation on the lowest energy 
states in the spectrum from the lowest of the large frequency {u)x) states is proportional 
to the square of the coupling strength divided by the energy difference by second order 
perturbation theory. This tells us that the one-dimensional limit is approached in the 
bottom of the spectrum with increasing deformation. This low-energy spectrum is 
much simpler and much less dense than that of the two-dimension spherical oscillator. 
Note that the spectrum is still two-fold degenerate due to the time-reversal symmetry, 
but now it is practically back to the standard equidistant scheme typical of harmonic 
confinement. 

4. Discussion 

We have shown that the single-particle spectrum for particles that are subject to a 
Rashba spin-orbit coupling in a deformed harmonic trap is very dependent on the 
parameters and that spectral gaps can be tuned by varying the Rashba coupling 
strength, a perpendicular Zeeman field and by changing the deformation of the 
trapping potential in the plane where the Rashba coupling acts. In particular, we 
find that the effect of deformation can be very different depending on whether the 
ratio of the trapping frequencies in the plane is equal to ratios of small integers or 
if it is closer to an irrational number. This can change the spectrum from having 
an intricate structure with several characteristic gap sizes into a more evenly spaced 
distribution that approximates a continuum. 

In this study we have not taken interactions into account. However, it is certainly 
possible to infer some potential effects of our results on a system with more particles 
and with interactions. The simple case of a non-interacting Fermi gas becomes very 
intriguing with the spectra presented here since a Fermi energy that is located in one 
of the pockets would be more stable than if it occurs in a region of higher state density. 
With the possibility of tuning this by external means we can explore this physics in a 
very direct manner. If we consider a weak interaction in a mean-field picture like the 
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BCS approach, the Rashba coupled case is extremely interesting as has been discussed 
recently [501 HD 1121 [21IM1 US . From BCS theory in finite systems and/or in a trap 
we know that one needs a critical coupling strength to obtain a superfluid state (in 
free space the critical value goes to zero). This critical value depends intimately on 
the density of states and the tunability demonstrated here shows that one can probe 
exotic superfluid states by using the deformation. 

For neutral single-species ultracold atoms the interactions are of van der Waals 
type and for the low experimental densities this can be considered of essentially zero- 
range. For a system of fermions with two internal hyperfine states this implies that 
one only has interaction in the spin singlet channel. The results presented here can 
be generalized to the two-particle case and one may then ask for the probability to 
find two fermions in the singlet state as the Rashba coupling strength, the Zeeman 
field and the deformation is varied. This will then tell us how the interactions are 
suppressed by the trapping potential in the presence of a spin-orbit coupling. In fact 
this is the first step toward the study of few-body physics in trapped system with 
non-abelian gauge potentials, and this is a direction that we will explore in the near 
future. Some open questions are the spectrum of two and three fermions in a trap 
with short-range interactions and its dependence on the deformation of the external 
confinement. Recent work has explored the case of few-boson problems in gauge 
fields in order to study strongly correlated dynamics akin to quantum Hall systems 
[12 133] . We expect that fermions in different trapping geometries should also give rise 
to interesting strongly correlated states when subjected to external gauge potentials. 
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